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ABSTRACT 

Context. Strontium has proven itself to be one of the most important neutron-capture elements in the study of metal-poor stars. Thanks 
to the strong absorption lines of Sr, they can be detected even in the most metal-poor stars and also in low-resolution spectra. However, 
we still cannot explain the large star-to-star abundance scatter we derive for metal-poor stars. 

Aims. Here we compare Galactic chemical evolution (GCE) predictions with improved abundances for Sr i and Sr n, including updated 
atomic data, to evaluate possible explanations for the large star-to-star scatter at low metallicities. 

Methods. We have derived abundances under both local thermodynamic equilibrium (LTE) and non-LTE (NLTE) for stars spanning 
a large interval of metallicities, as well as a broad range of other stellar parameters. Gravities and metallicities are also determined in 
NLTE. We employed MARCS stellar atmospheres and MOOG for the LTE spectrum synthesis, while MAFAGS and DETAIL were 
used to derive the NLTE abundances. We verified the consistency of the two methods in LTE. 

Results. We confirm that the ionisation equilibrium between Sri and SrII is satisfied under NLTE but not LTE, where the difference 
between neutral and ionised Sr is on average ~ 0.3 dex. We show that the NLTE corrections are of increasing importance as the 
metallicity decreases. For the stars with [Fe/H] > -3, the Sri NLTE correction is ~ 0.35/0.55 dex in dwarfs/giants, while the Srn 
NLTE correction is < ±0.05 dex. 

Conclusions. On the basis of the large NLTE corrections to Sr i, Sr i should not be applied as a chemical tracer under LTE, while it 

is a good tracer under NLTE. Sr n, on the other hand, is a good tracer under both LTE and NLTE (down to [Fe/H] 3), and LTE 

is a safe assumption for this majority species (if the NLTE corrections are not available). However, the Sr abundance from Srn lines 
depends on determining an accurate surface gravity, which can be obtained from the NLTE spectroscopy of Fe lines or from parallax 
measurements. We could not explain the star-to-star scatter (which remains under both LTE and NLTE) by the use of the Galactic 
chemical evolution model, since Sr yields to date have been too uncertain to draw firm conclusions. At least two nucleosynthetic 
production sites seem necessary to account for this large scatter. 
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1. Introduction 

Strontium is one of the two neutron-capture elements (namely 
strontium and barium - Ba) that intrinsically show very strong 
absorption lines even in metal-poor stars. At solar metallicity 
Sr is synthesised by a variety of nucleosynthetic sources includ- 
ing the weak slow neutron-capture (s-) process that occurs in 
massive stars (e.g. Heil et al. 2009; Pignatari et al. 2010) and in 
AGB stars (Travaglio et al. 2004). In comparison, the produc- 
tion of Ba is dominated by the s-process occurring in low-mass 
AGB stars (Kappeler et al. 1989; Busso et al. 1999; Sneden et al. 
2008). This picture changes at low metallicity, where Ba may 
be formed by a main rapid neutron-capture process and Sr by 
a charged particle process (Hoffman et al. 1997). Additionally, 
very metal-poor rapidly rotating massive stars might be signif- 
icant producers of Sr and Ba via the s-process (Pignatari et al. 



2008; Chiappini et al. 201 1; Frischknecht et al. 2012). This may 
modify the presumption of pure r-process patterns in ultra metal- 
poor (UMP; -5 < [Fe/H] < -4 cf. Beers & Christlieb 2005) 
stars. For this reason, disentangling the nucleosynthetic origin 
of Sr and Ba in metal-poor stars would help us understand the 
formation and evolution of the early Galaxy. 

Only the 4077 A Sr II line remains detectable in both dwarfs 
and giants both in high- and low-resolution spectra of metal-poor 
stars. Studying this line thus provides unique insight into the be- 
haviour of neutron-capture elements at all metallicities and spec- 
tral resolutions, ranging from the low-resolution LAMOST sur- 
vey to the high-resolution Gaia-ESO survey 1 . Clearly, accurate 
abundances are needed to fully comprehend the chemical evo- 

1 The Srn line is detectable if a blue setting is used for Gaia-ESO 
follow-up observations. 
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lution of Sr. This means that the effects of non-local thermody- 
namic equilibrium (NLTE) and deviations from hydrostatic equi- 
librium (3D) must be taken into account in element abundance 
calculations, if we want to extract the correct information from 
the future surveys' large flow of data. Clearly, such calculations 
are a challenge, as was recently demonstrated for O (Gonzalez 
Hernandez et al. 2010), Ca (Spite et al. 2012), Fe (Bergemann 
et al. 2012a), and for Ba (Dobrovolskas et al. 2012), while 
Bonifacio et al. (2009) provide 3D corrections for a large number 
of elements for dwarfs. The estimates of 3D effects for the Srn 
resonance lines have been provided by Collet et al. (2007), who 
performed LTE calculations with 3D radiative-hydrodynamics 
simulations of stellar convection for the metal -poor stars . The 3D 
LTE corrections are of the order of -0.15 dex, with respect to ID 
LTE. The NLTE abundances of Sr were reported in a few studies, 
such as those of Belyakova & Mashonkina (1997), Mashonkina 
et al. (1999), Andrievsky et al. (2011), and Bergemann et al. 
(2012a). For the Srn resonance line, the NLTE abundance cor- 
rections are not large, typically within +0.2 dex 2 , however, they 
are sensitive to variations in the stellar parameters (Bergemann 
et al. 2012a). Full 3D NLTE calculations for Sr have not been 
performed yet. 

Using the NLTE technique presented in Bergemann et al. 
(2012a), we have now expanded the stellar sample in order to 
study the chemical evolution of Sr in the Galaxy. We have de- 
rived NLTE Sr abundances and NLTE stellar parameters for a 
sample of 21 stars, plus comparison samples (51 very metal- 
poor stars from Francois et al. 2007 and Bonifacio et al. 2009). 
We also include the predictions of the homogeneous chemi- 
cal evolution model for the Galactic halo of Chiappini et al. 
(2008), computed with the most up-to-date Sr yields available 
in the literature. Given the still large uncertainties on the stellar 
yields, the goal of comparing the data and the chemical evolu- 
tion model is to give a first impression of how far the available 
stellar yields are from explaining the data. A comparison with in- 
homogeneous chemical evolution models, more suitable to low- 
metallicity environments, is beyond the scope of the present pa- 
per, see Cescutti (2008); Cescutti & Chiappini (2010). The paper 
is structured as follows. Sections 2 and 3 describe the observa- 
tions, stellar parameters, and NLTE calculations, respectively. 
Sections 4 and 5 present the results, yields, and the chemical 
evolution model. Conclusions can be found in Sect. 6. 

2. Sample and data reduction 

A sub-sample of stars (marked 'u' in Sect. 4) was taken 
from Hansen et al. (2012). These stars were observed with 
UVES/VLT (Dekker et al. 2000, R > 40000) between 2000 and 
2002, and their spectra were reduced with the UVES pipeline 
(v. 4.3.0). The spectra have a signal-to-noise ratio, S/N, > 100 
at 3200 A. Two stars (see 'h' in the following) were observed 
with HIRES/Keck (R ~ 60000, Vogt et al. 1994). Their spectra 
are of similar quality as the UVES spectra, and they have been 
retrieved from the HIRES reduced data archive. The pipeline- 
reduced data were wavelength-shifted, co-added, and had their 
continua normalised before our analysis. For further details we 
refer to Hansen et al. (2012). 

Three stars (HD 134169, HD 148816, HD 184448) were ob- 
served with the FOCES echelle spectrograph at the 2.2 m tele- 
scope of the CAHA observatory on Calar Alto, during 1999 and 

2 For very metal-poor (VMP; -3 < [Fe/H] < -2) stars the corrections 
are less than about +0.08 and slightly greater for extremely metal-poor 
stars ([Fe/H] < -3). 



2000, and were kindly made available to us by T. Gehren. The 
spectra have a resolution of ~ 60000 and an S /N of ~ 200 near 
5000 A. More details on the observations and data reduction can 
be found in Gehren et al. (2004, 2006). 

We selected the stars according to the following criteria: 
1) the observations cover the spectral range of the Sri 4607 A 
line, 2) accurate photometry is available, and 3) the stars cover 
a broad stellar parameter space to test the Sri and Srn abun- 
dance behaviour at different temperatures, gravities, and metal- 
licities. Our sample thus consists of 21 dwarf, sub-giant, and gi- 
ant stars. Here we have disregarded carbon enhanced metal-poor 
stars (CEMP-s/CEMP-no — with and without s-process over- 
abundances). Since the CEMP-s stars tend to have very large s- 
process abundances, we are slightly biased against high Sr abun- 
dances at low metallicity. 

For comparison, we include extremely metal-poor stars from 
Francois et al. (2007) and Bonifacio et al. (2009). For details 
about the observed data, we refer to these publications. 

3. Methods 

In this work, NLTE effects are accounted for in the determination 
of basic stellar parameters (surface gravity and metallicity), as 
well as Sr abundances. We describe the analysis in detail below. 

3.1. Model atmospheres 

All calculations in this work were performed with classical ID 
LTE plane-parallel model atmospheres. We used the MAFAGS 
(Grupp 2004a,b) and MARCS (Gustafsson et al. 2008) mod- 
els, which are both well-adapted to analysis of late-type stars. 
The model atmosphere codes adopt slightly different prescrip- 
tions for the convective flux transport and background opacity. 
However, comparison of the model T(r) relation showed that 
these differences are very small, and reveal themselves only in 
optically thick layers, where the treatment of convection is im- 
portant (see Fig. 1). A comparison of the MAFAGS and MARCS 
models was presented in Bergemann et al. (2012b), where we 
showed that the differences in Fe abundances obtained with dif- 
ferent model atmosphere codes are very small, typically within 
0.05 dex. We comment on this further in Sect. 3.3. 

3.2. Stellar parameters 

Stellar parameters for the selected sample of stars were taken 
from (Hansen et al. 2012, 'H') and (Bergemann & Gehren 2008; 
Bergemann et al. 2012b, - 'BG';'B'), giving preference to pa- 
rameters determined with IR photometry and parallaxes. The 
'H', 'BG', and 'B' indicate the source of the temperature here- 
after. A brief description of these data follows. 

The effective temperatures were derived from several differ- 
ent colour indices and calibration methods (Alonso et al. 1996, 
1999; Ramirez & Melendez 2005; Masana et al. 2006; Onehag 
et al. 2009; Casagrande et al. 2010). We chose temperature cal- 
ibrations that fall in the middle of the probed calibration ranges 
(see Hansen et al. 2012, for details). The reddening values were 
taken from Schlegel et al. (1998), and since all the reddening 
values are much lower than 0. 1 mag (Bonifacio et al. 2000), we 
have not applied their corrections to these values. Our chosen 
stellar parameters are also consistent with the effective temper- 
atures (r e ff) determined from the ID fitting of Balmer profiles 
by Gehren et al. (2004, 2006). For the few stars (HD 19445, 
HD 142038, G 64-12) we have in common with Gehren et al. 
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to +0.3 dex. For the extremely metal-poor (EMP; -4 < [Fe/H] 
< -3) giant stars, we obtain a maximum NLTE correction of 
Anlte [Fe/H] of +0.3 dex. 
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Fig. 2: Comparison of LTE and NLTE surface gravities (top) and 
metallicities (bottom). 



Fig. 1: Comparison of MARCS (dotted) and MAFAGS (solid) 
model atmospheres for selected stellar parameters at two differ- 
ent metallicities ([Fe/H] = -1.2, -2.4). 



(2004, 2006), they derive T eS = 5985, 5773, and 6407 K, re- 
spectively, with an rms offset of 10 K from our values. This con- 
firms the agreement between the B aimer r e ff scale from Gehren 
et al. (2004, 2006) and the method described in Alonso et al. 
(1996) that we used here. Balmer line T e s's (ID) from the same 
reference agree with the values we adopted for HD 134169, HD 
148816, and HD 184448. The parameters for HD 122563 and G 
64-12 are those from Bergemann et al. (2012b). 

For the stars with parallax measurements, the surface gravity 
was calculated using the classical formula that relates mass, tem- 
perature, magnitude and parallax to gravity. Masses and bolo- 
metric corrections were taken from Nissen et al. (1997, 2002, 
2007). Metallicities were then initially estimated in LTE, and 
the effects of NLTE were taken into account by applying NLTE 
abundance corrections to Fe. These corrections were computed 
for the adopted Fe i line list by interpolation in the Fe NLTE 
grid presented by Lind et al. (2012). We note that for stars with 
metallicity [Fe/H] > -2, the systematic difference between LTE 
abundances of Fe i and Fe n is not large (< 0. 1 dex). The Sr NLTE 
corrections also stay within 0.05 - 0.07 dex (see Sect. 4) consis- 
tent with the results in Bergemann et al. (2012a). The effect of 
NLTE becomes very important for Fe in very metal-poor stars 
with [Fe/H] < -2.5. In particular, most of the stars from our 
comparison samples (Francois et al. 2007; Bonifacio et al. 2009) 
are subject to NLTE metallicity corrections ranging from +0.2 



For the other stars, which do not have known parallaxes, the 
initial estimate of log g and [Fe/H] was obtained from the LTE 
ionisation equilibrium of Fe, and Fe NLTE corrections were ap- 
plied to both gravity and metallicity. To compute the NLTE cor- 
rection for surface gravity, we used the approximate Anlte log g 
- ANLTE[Fe/H] calculations from Lind et al. (2012, their Sect. 
3). The changes in \ogg due to NLTE effects are significant and 
have a clear effect on the [Sr/Fe] ratios derived from the gravity- 
sensitive Sr n lines, compared to the changes in metallicity. The 
NLTE gravity corrections reach up to ~ +0.8 dex for the most 
metal-poor giants in our sample. The influence of T e g on Srn is 
minor. By changing the temperature with its uncertainty, the Sr n 
abundance changes by 0.01-0.05 dex, while the change in Sri is 
much larger (see Sect. 4.1). The NLTE and LTE stellar parame- 
ters are compared in Fig. 2. 

As seen from Fig. 3 the temperatures for seven stars taken 
from three different sources (Casagrande et al. 2010; Bergemann 
et al. 2012a; Hansen et al. 2012) agree within 40 K in most cases, 
and for a few stars (e.g. HD106038) the difference between the 
two IRFM methods (Alonso et al. 1996; Casagrande et al. 2010) 
is 170 K. This difference is within the combined errors, if we in- 
clude systematic errors, as well as the uncertainty on E(B - V). 
The dwarfs from the comparison sample (Bonifacio et al. 2009) 
had their temperatures estimated from Ho- line profile fitting, 
which yielded values in good agreement with those determined 
from Alonso et al. (1996) calibrations. We refer the reader to 
Bonifacio et al. (2009) and Sbordone et al. (2010) for further de- 
tails on Ha line profile fitting. The giants from the second com- 
parison sample (Francois et al. 2007) had their temperatures de- 
termined by the use of broad range photometry calibrations from 
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Fig. 3: Comparison of effective temperatures determined with 
different methods for seven stars. The legend indicates the 
original paper (Casagrande et al. 2010; Hansen et al. 2012; 
Bergernann et al. 2012a, - yellow '*', blue diamond, and green 
square, respectively) from which the temperatures have been 
taken. 



Alonso et al. (1999). This is the same method we applied here, 
and the differences/offsets between our temperature determina- 
tions and those made for the comparison samples are minimal. 

The microturbulence velocity was fixed by requiring that Fe i 
lines yield the same abundances regardless of their equivalent 
width. In our parameter space, the NLTE corrections to micro - 
turbulence (Lind et al. 2012) are smaller than the formal uncer- 
tainty of (+0.15 km/s) and were not considered here. 

The main error in our photometric temperatures comes from 
reddening (< +0.05 mag). For Balmer lines, the temperature 
error is largely internal and is determined from the profile fit- 
ting. This applies to the six stars marked by a 'B' in the fol- 
lowing. The error in the adopted \ogg values is dominated by 
that of parallaxes (< +1.0"). The errors in [Fe/H] and mi- 
croturbulence are assumed to be 0.1(5)dex and 0.15 km/s, re- 
spectively. After propagating all stellar parameter uncertain- 
ties, we adopted a common set of uncertainties for our stars 
of (r eff /logg/[Fe/H]/£): ±100K/0.2dex/0.1 dex/0.15km/s. For 
the most metal-poor stars, we find slightly higher values: 
±100 K/0.25 dex/0. 15 dex/0. 15 km/s. These errors are internal to 
our method. The differences between LTE and NLTE stellar pa- 
rameters highlights that systematic errors in LTE are larger than 
the internal errors. 

3.3. Sr abundance determinations 

The NLTE statistical equilibrium calculations for Sr were per- 
formed with the revised version of the DETAIL code (Butler 
& Giddings 1985). The new model atom of Sr and other re- 
lated aspects of the NLTE calculations are described in detail in 
Bergernann et al. (2012a). The NLTE effects on the Sri lines are 
primarily caused by over-ionisation, which leads to systemati- 
cally higher NLTE abundances compared with LTE, especially 
for more metal-poor and hotter stars. In contrast, deviations from 
LTE in the Sr n lines are largely driven by strong line scattering. 
As a consequence, the differences between the LTE and NLTE 
abundances may be positive or negative, depending on tempera- 
ture, gravity, and metallicity of a star. 
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Fig. 4: Effect of elastic H I collisions on the 4077 A lines ex- 
pressed through different C6 values (red dashed and black solid 
line) - compared to observations (black dots) of a dwarf (HD 
106038) and a metal-poor giant (HD 122563) star. 



The LTE and NLTE abundances of Sr in the selected metal- 
poor stars were determined as follows. The LTE abundances are 
synthesised with MOOG, while NLTE synthesised abundances 
are derived using the SIU code (Reetz 1999), where the NLTE 
departure coefficients are computed with DETAIL. In both cases 
we apply the same atomic data. We note that the LTE abundances 
from SIU and MOOG agree within 0. 1 dex. This difference is 
a combination of local continuum placement (which can be up 
to 0.05 dex when set locally by eye) and of the different syn- 
thesis codes, since the model atmospheres are very similar, as 
seen from Fig. 1, they almost do not contribute to this differ- 
ence. Only very deep in the atmosphere (close to r = 1) do the 
MARCS and MAFAGS models differ due to a more efficient 
transport of convective flux in the MARCS models, which leads 
to a cooler atmosphere than in MAFAGS. For 18 stars, with a 
four-star overlap with the SIU analysis, we first determine LTE 
abundances by spectrum synthesis with the MOOG code. To ob- 
tain NLTE abundances, we then apply the NLTE corrections cal- 
culated with DETAIL and MAFAGS. We do not perform differ- 
ential abundance analysis with respect to the Sun. 

In particular, to remain consistent with our previous analy- 
sis, we use the van der Waals broadening and ^/-values from 
Bergernann et al. (2012a) (see Table 1). According to this study 
the damping constants for the Srn lines could be somewhat un- 
certain. Variation in log C(, by +0.35 (~ 20%) leads to a change 
in abundances by +0.15 dex for the 4077 A Sr n line in the metal- 
poor stars with [Fe/H]> -1.5 (see Fig. 4), but only has a minor 
effect at very low metallicity (~ +0.05 dex). The accuracy of 
the gf-value for the 4607 A Sr i line was critically evaluated by 
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Table 1 : Atomic data for Sr I and II 



HD 122563, Sr II 



Element 


A 


X 


l°gg/ 


log(y/N„) 


log C 6 




A 


eV 




rad cm 3 s _1 


cm 6 s 


Sri 


4607.33 


0.0 


0.283 


-7.53 


-31.2 


SrII 


4077.71 


0.0 


0.158* 


-7.81 


-32.0 



Notes. ( *' Total log gf value. Further details and hfs splitting can be 
found in Bergemann et al. (2012a). 
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Fig. 5: Synthetic spectra (solid lines) of the 4607 A Sr i line plot- 
ted on top of the observed spectrum of HD 106038 (dots). 



the NIST database. The uncertainties are less than one percent 
and yield accurate abundances even for metal-rich stars (see Fig. 
5), whereas the Sr n resonance lines at 4077 and 4215 A are too 
blended and strong to give any reliable information about the 
solar Sr abundance (see Fig. 6). We derive our Sr i abundances 
from the 4607 A line and the Sr n abundances from the 4077 A 
line, since all our stars have sub-solar metallicities. 

We tested the Abfind package in MOOG, which uses 
the measured equivalent widths (EW) to compute abundance 
by the curve-of-growth method. We find, however, that this 
method yields abundances in slight disagreement (+0.04 < 
[Sr/Fe]Ew-syn < ±0.25) with the results obtained by the full 
profile fitting using the Synth (synthesis) package of the same 
code (see the online material). Depending on the line properties, 
in particular equivalent width, the abundance is either over- or 
under-estimated (Table 2). This expected; the Sri line is very 
weak in metal-poor stars, and both Sr n lines are strong, sensi- 
tive to damping, and affected by blends. For example, the EWs 
determined by fitting Voigt and Gaussian profiles in IRAF ('/' 
separated entries in the 'EW' columns in Table 2) differ gener- 
ally by 5-10% (in a few cases, like the VMP dwarf HD 106038, 
the difference is approximately a factor of three larger). As a re- 
sult, abundances derived using the EWs may be discrepant by up 
to 0.4 dex (Fig. 7). However, except for one case (the EMP sub- 
giant HD 140283), the largest difference between EW and syn- 
thesis determined abundances are of the order of +0.1 dex (see 
online material). This value might be slightly over-estimated due 
to local continuum placement. 

Blends influence abundance determinations in both the 
metal-poor and metal-rich parts of our sample. Features in the 
blue and red wings of the 4077 Srn line, ~ 0.3 A away from the 
line centre, only vanish in the spectra of most metal-poor stars, 
[Fe/H] < -2. These blends are due to lanthanum, chromium, and 
dysprosium. The Srn line at 4215 A is blended by the two Fei 




Fig. 6: Comparison of synthetic profiles (lines) of the Srn line 
to observations of HD 122563 (black dots). The syntheses have 
been computed with (black solid line) and without (red dashed 
line) the Fe i blends for the 4215 A Srn line. 
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Fig. 7: Difference in Sr abundances between EW and synthesis. 
Dwarfs/giants are shown as small red/big blue filled circles, re- 
spectively. 



lines at 4215.42 and 4216.18 A. The former is rather strong and 
clearly distorts the shape of the Sr n profile, as directly seen in 
the very high-resolution spectra. Figure 6 shows that even for 
the very metal-poor giant HD 122563 with [Fe/H]NLTE = -2.5, 
the abundance is over-estimated by 0.14 dex (when synthesised) 
if this blend is not taken into account. Therefore, the abundance 
determined from the 4215 A Srn line is subject to a systematic 
uncertainty, and the blend is sensitive to T e g , log g, and metallic - 

ity. 

On these grounds, we do not include the 4215 A Srn line 
in the abundance calculations. Furthermore, for the comparison 
samples from Francois et al. (2007) and Bonifacio et al. (2009), 
only measurements for the 4077 A Sr n line are available. To be 
consistent and to retain the full sample size, we therefore only 
use the 4077 Srn line. 

The Sr i line is generally weak, and regardless of the profile 
fitted to this line, the EW-converted abundances tend to be larger 
than the synthesised abundances. This might be due to an iron 
blend in the red wing of this neutral strontium line (see Fig. 5, 
where a stronger Fe line blends into the red wing of the Sr i line). 

The LTE abundances derived for our sample were calculated 
with LTE stellar parameters (as described in Sect. 3.2) using the 
ID LTE synthetic spectrum code MOOG to synthesise spectra 
for these stars. The EWs for the dwarf comparison sample were 
taken from Bonifacio et al. (2009). We measured the EWs for the 
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Table 2: Example of LTE [Sr/Fe] abundances for a dwarf, sub-giant, and giant from the measured equivalent widths (EW) and from 
synthesis (synt). The measured Gauss/Voigt EWs (G/V) are given in mA , together with the stellar parameters. 



Star 


4077 EW 


EWG/V 


4077 synt 


4607 EW 


EWG/V 


4607 synt 


T c ff 




[Fe/H] Lr£ 








[mA] 






[mA] 




[K] 






[km/s] 


HD 10603 8 


0.28/0.65 


186.1/271.0 


0.5 


0.43/0.47 


13.1/14.1 


0.33 


5950 


4.33 


-1.48 


1.1 


HD140283 


-0.44/ - 0.33 


75.8/79.8 


-0.15 








5777 


3.70 


-2.58 


1.5 


HD 122563 


-0.34/ - 0.09* 


158.8/184.0 


-0.05 


-0.54* 


2.9 


< -0.6 


4665 


1.65 


-2.50 


1.8 



Notes. w Value uncertain 



remaining stars. The LTE Sr II abundances for the comparison 
samples (dwarfs and giants) were redetermined with MARCS 
models and MOOG, using the solar abundance, from Anders 
& Grevesse (1989), which we have adopted for this study. The 
LTE abundances we calculated agree within 0.05 -0.1 dex with 
those published in (Francois et al. 2007, F07) and (Bonifacio 
et al. 2009, B09), and only for a handful of stars is the difference 
greater than 0.2 dex. 

The correction for NLTE effects is most conveniently per- 
formed by differentiating LTE and NLTE curves-of-growth at 
a given line strength. The method we adopted thus relies on 
the determination of EWs and subsequent translation to LTE 
and NLTE abundances. From the comparison sample (the "First 
Stars" samples), we only have these EWs. When determin- 
ing EWs by profile fitting with a single component, unresolved 
blends may play a role, so we derived both LTE and NLTE 
EW-based abundances and compared them to synthesised abun- 
dances to assess the impact the single-line assumption have on 
the final abundances. For all the tests related to atomic data, 
blends and profile fitting, we have maintained one set of stellar 
parameters for each star (those listed in Table 3). Thus, the re- 
sulting difference in abundance is an expression of uncertainties 
in the atomic data, unknown blends, and continuum placement. 

From Fig. 7 we estimate that the Sr abundances from the 
metal-rich stars might be overestimated when using the EW 
method, while Sr in the metal-poor stars will be overestimated 
with < 0.32 dex. Another part of the over-/under-estimation can 
be assigned to the line profile fitting and continuum placement. 
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Fig. 8: Sr LTE (open circles) and NLTE abundance differences 
(filled big blue/small red circles for giants/dwarfs) as a function 
of [Fe/H] (top panel) and logg (bottom panel). 



4. Results: LTE vs NLTE 

The results for our sample are summarised in Table 3, which pro- 
vides the mean of synthesised abundances for Sr i and Sr n. Their 
differences are shown in Fig. 8 as a function of stellar [Fe/H] and 
log g. The LTE and NLTE line by line abundances from both EW 
and synthesis are also given in the online Table A. 1 (for our sam- 
pie). 

The LTE approximation fails to establish the ionisation bal- 
ance of Sri and Srn. Figure 8 shows that the offset between the 
two ionisation stages is about 0.2 dex for dwarfs, but it increases 
up to 0.5 dex for giants. This difference is mainly caused by 
the progressively increasing systematic error in the LTE abun- 
dance inferred from the Sr i line, which shows NLTE abundance 
corrections of up to 0.5 dex at low metallicity and low gravity 
(online Table A. 1). Although the NLTE effects on the resonance 
Srn lines are not significantly pronounced, they depend on stel- 
lar parameters, particularly on the [Fe/H] or, equivalently, on the 
Sr abundance itself (see discussion in Sect. 3.3). In Table 3, we 
see that the LTE abundances obtained from the Sr n lines can be 
over- or under-estimated by up to 0.1 dex, which may introduce 
a spurious systematic trend or, more likely give rise to a larger 



line-to-line scatter. If the gravities are derived from parallaxes, 
their values will be the same in both LTE and NLTE; however, 
when derived from ionisation equilibrium, the NLTE corrected 
log g CO r r will differ from the LTE Fe i based log g. For our stel- 
lar sample, even though the line-to-line scatter is clearly smaller 
under NLTE, the star-to-star scatter almost remains the same un- 
der both LTE and NLTE. The results for the comparison samples 
(the 'First Stars ' samples - F07, and B09) are shown in Table 4. 

Figure 9 shows the [Sr/Fe] ratios as a function of [Fe/H]. In 
the figure, the error bars are the total, propagated uncertainties, 
computed as described below. The exceptions are HD 3567, HD 
19445, HD 122563, and HD 126587 with upper limits on the 
abundance from the 4607 A Sr i line. 

4.1. Uncertainties 

A number of test calculations varying the input parameters in 
the spectrum synthesis were performed for the two representa- 
tive stars with the same metallicity: HD 106038 (dwarf) and HD 
74462 (giant). In particular, we are interested in the sensitivity of 
the abundances to the model atmosphere parameters: T e s, \ogg, 
[Fe/H], and microturbulence (Table 5). As seen from this table, 
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Table 3: Stellar parameters and derived Sr abundances ± standard deviation for the selected sample of stars. Strontium abundances 
with superscript (1) are derived only from the 4077 A Srn line. Parameters with subscript have been corrected for NLTE effects. 
Stars with an V superscript have had their abundances derived in SIU. 



Star 


n ±<r 




log s 


log gcorr. 


[Fe/Hl/TF 


fFe/Hljv/TF 

L 1 J/VLi 




[Sr/Fel/ t£ 


rSr/FelAz/TF 


Comment 




[mas] 


[K] 










km/s 








HD 3567 


9.57+1.38 


6035 


4.08 


4.08 


-1.33 


-1.29 


1.5 


-0.03 ± 0.8* 


0.06 ± 0.8* 


[u,H] 


HD 19445 


25.85+1.14 


5982 


4.38 


4.38 


-2.13 


-2.10 


1.4 


0.13 ±0.8* 


0.16 + 0.8* 


[u,H] 


HD 106038 


9.16+1.50 


5950 


4.33 


4.33 


-1.48 


-1.45 


1.1 


0.42+0.12 


0.45 ± 0.21 


[u,H] 


HD 121004 


16.73+1.35 


5711 


4.46 


4.46 


-0.73 


-0.71 


0.7 


0.18 + 0.04 


0.26 + 0.17 


[u,H] 


HD 122196 


9.77+1.32 


6048 


3.89 


3.89 


-1.81 


-1.75 


1.2 


0.24 (1) 


0.19 m 


[u,H] 


HD 134169 


16.80+1.11 


5930 


3.98 


3.98 





-0.86 


1.8 


-0.05 + 0.14 


-0.06 + 0.13 


[f,s,BG] 


HD 140283 


17.16+0.68 


5777 


3.70 


3.70 


-2.58 


-2.38 


1.5 


-0.15 (1) 


-0.37 (1) 


[u,H,B] 


HD 148816 


24.34+0.90 


5880 


4.07 


4.07 




-0.78 


1.2 


-0.13 + 0.18 


-0.13 + 0.17 


[f,s,BG] 


HD 184448 


19.16+0.63 


5765 


4.16 


4.16 





-0.43 


1.2 


0.00 ± 0.21 


-0.01+0.21 


[f,s,BG] 


G 64-12 


0.57+2.83 


6464 


4.30 


4.30 


-3.24 


-3.12 


1.5 


0.00 (1) 


0.17 (1) 


[u,B] 


G 64-37 


2.88+3.10 


6494 


3.82* 


4.23 


-3.17 


-3.00 


1.4 


0.08 (1) 


0.17« 


[u,H] 


HD 122563 


4.22+0.35 


4665 


1.65 


1.65 


-2.60 


-2.50 


1.8 


-0.23 ± 0.8* 


-0.12 + 0.8* 


[u,B] 


HD 175305 


6.18+0.56 


5100 


2.70 


2.70 


-1.38 


-1.34 


1.2 


-0.12 + 0.32 


0.04 + 0.1 


[h,H] 


BD -133442 




6450 


4.20 


4.42 


-2.56 


-2.47 


1.5 


0.30 (1) 


0.21* 1 ' 


[u,H] 


CS 30312-059 




5021 


1.90 


2.41 


-3.06 


-2.89 


1.5 


0.50 (1) 


0.31<" 


[h,H] 


CS 31082-001 




4925 


1.51 


2.05 


-2.81 


-2.63 


1.4 


0.70 (1) 


0.60<" 


[h,H] 


HD 74462 




4590 


1.84 


1.98 


-1.48 


-1.43 


1.1 


-0.25 ± 0.35 


-0.14 + 0.03 


[h,H] 


HD 126238 




4900 


1.80 


2.02 


-1.92 


-1.85 


1.5 


-0.17 + 0.24 


0.01 ± 0.06 


[u,H] 


HD 126587 




4950 


1.90 


2.36 


-3.01 


-2.86 


1.65 


0.23 ± 0.8* 


0.22 ± 0.8* 


[u,H] 


HE 0315+0000 




5050 


2.05 


2.47 


-2.81 


-2.67 


1.7 


0.39 (1) 


0.23 (1) 


[u,H] 


HE 1219-0312 




5100 


2.05 


2.58 


-2.99 


-2.81 


1.65 


0.29 (1) 


0.12<» 


[u,H] 



Notes. w Value from ionisation equilibrium 

(b,bg,h) x em p era t ur e from (Bergemann et al. 2012b, B), (Bergemann & Gehren 2008, BG), and (Hansen et al. 2012, H), respectively. 
(«,/,/.) Observed spectra from: " UVES/VLT, / FOCES/Calar-Alto Gehren et al. (2004, 2006), * HIRES/Keck 

(*) Weighted average, which includes an upper limit that was given half weight. The large uncertainty is a reciprocal square root of the summed 
weights. 



the Srn line at 4077 A is very strong and mainly sensitive to 
log g and In comparison, the neutral Sr abundance from the 
weaker 4607 A line is sensitive to temperature and [Fe/H], but 
almost not affected by log g and £. 

Several of the Sr i abundances presented in the online Table 
A.l are only upper limits, which is why we applied a weighted 
average and a corresponding reciprocal square root of the 
summed weights (see Taylor 1997). The final [Sr/Fe] abun- 
dances and their uncertainties are shown in Table 3. The upper 
limits of Sr i were given half the weight of Sr n. This overesti- 
mates the uncertainty a bit. However, with only two measure- 
ments, this approach seems sensible. In some stars only Srn 
could be measured, so the abundance is based on only one trust- 
worthy line. To estimate the uncertainty in this case we made 
independent measurements of the 4077 Sr n line, and found the 
derived abundances to be consistent to within 0.0 - 0.1 dex. An 
average value of +0.05 dex was adopted instead of the standard 
deviation otherwise applied to the stars with two detectable lines. 

We summarise that the assumption of LTE especially in 
metal-poor, low-gravity stars will, in addition to not fulfilling 
Sr ionisation balance, also introduce a weak spurious trend of 
Sr abundances with metallicity. The Sri resonance line would 
consistently under-estimate the Sr abundance trend, whereas the 
abundances obtained from the subordinate Sr n lines would be 
systematically too large by ~ 0.05 dex. Whereas the abundances 
derived from the Sri line are significantly affected by NLTE 
line formation, Srn lines are less so. However, NLTE effects 
for Fe must be accounted for in the spectroscopic gravity deter- 



Table 5: Uncertainties in the individual LTE Sr I and II abun- 
dances in a giant (HD 74462) and a dwarf (HD 106038) star. 



HD 74462: [Sr/Fe] 


0.1 


-0.5 


Parameter/line [A] 


4077 


4607 


T ±100 


±0.03 


±0.17 


logg ±0.2 


±0.04 


±0.01 


[Fe/H] ±0.1 


±0.12 


±0.09 


£ ±0.15 


±0.07 


±0.01 


Propagated uncertainty 


±0.15 


±0.19 


HD 106038: [Sr/Fe] 


0.65 


0.33 


T±100 


±0.04 


±0.09 


logg ±0.2 


±0.17 


±0.01 


[Fe/H] ±0.1 


±0.22 


±0.09 


£±0.15 


±0.13 


±0.01 


Propagated uncertainty 


±0.3 


±0.13 



mination in order to derive accurate abundances for Sr II lines. 
Furthermore, we note that the large star-to-star scatter found in 
LTE abundance studies remains under NLTE even at extremely 
low metallicities. 

Our analysis of the chemical evolution of Sr in NLTE is dif- 
ferent from previous studies. Andrievsky et al. (2009) has al- 
ready performed NLTE calculations for Sr; however, their stel- 
lar parameters were determined assuming LTE. These have a 
measurable impact on the Sr abundances. The difference in 
the overall trend of [Sr/Fe] with metallicity (Bergemann et al. 
(2012a) and our Fig. 9) is detectable, and is best seen in Fig. 
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Table 4: Basic parameters and EW calculated Sr abundances from the 4077A line for the comparison stars. Top part contains dwarf 
stars from the comparison sample, while lower part shows giants from the second comparison sample. Details and subscripts are 
described in Table 3. 



Star 


T c ff 
[K] 


logg 


^corrected 


[FefH] LTE 


[Fe/H] NLTE 


[km/s] 


[Sr/Fe] ir£ 


[Si/Fe] NLTE 


Dwarf sample"' J * jy 


















BS 16023-046 


6364 


4.50 


4.69 


-2.97 


-2.90 


1.3 


-0.20 


-0.17 


BS 16076-006 


5199 


3.00 


3.29 


-3.81 


-3.70 


1.4 


0.67 


-0.64 


BS 16968-061 


6035 


3.75 


4.04 


-3.05 


-2.94 


1.5 


-1.59 


-1.58 


BS 17570-063 


6242 


4.75 


4.87 


-2.92 


-2.87 


0.5 


0.03 


-0.02 


CS22177-009 


6257 


4.50 


4.67 


-3.10 


-3.03 


1.2 


-0.15 


-0.12 


CS22888-031 


6151 


5.00 


5.09 


-3.30 


-3.26 


0.5 


0.05 


0.09 


CS22948-093 


6356 


4.25 


4.53 


-3.30 


-3.19 


1.2 


-0.08 


-0.01 


CS22953-037 


6364 


4.25 


4.48 


-2.89 


-2.80 


1.4 


-0.45 


-0.46 


CS22965-054 


6089 


3.75 


4.06 


-3.04 


-2.92 


1.4 


-2.09 


-2.17 


CS22966-011 


6204 


4.75 


4.87 


-3.07 


-3.02 


1.1 


0.88 


0.95 


CS29499-060 


6318 


4.00 


4.26 


-2.70 


-2.60 


1.5 


-0.63 


-0.70 


CS29506-007 


6273 


4.00 


4.27 


-2.91 


-2.80 


1.7 


-0.49 


-0.50 


CS29506-090 


6303 


4.25 


4.46 


-2.83 


-2.75 


1.4 


0.33 


0.27 


CS295 18-020 


6242 


4.50 


4.65 


-2.77 


-2.71 


1.7 


0.08 


0.05 


CS295 18-043 


6432 


4.25 


4.53 


-3.20 


-3.09 


1.3 


— 


— 


CS29527-015 


6242 


4.00 


4.34 


-3.55 


-3.41 


1.6 


0.12 


0.22 


CS30301-024 


6334 


4.00 


4.27 


-2.75 


-2.64 


1.6 


-0.41 


-0.47 


CS30339-069 


6242 


4.00 


4.28 


-3.08 


-2.97 


1.3 


0.43 


0.40 


CS31061-032 


6409 


4.25 


4.45 


-2.58 


-2.50 


1.4 


-0.48 


-0.54 


Giant sample'' Sftul 


















BD+17:3248 


5250 


1.40 


1.97 


-2.07 


-1.88 


1.5 


0.00 


-0.09 


BD-18:5550 


4750 


1.40 


1.95 


-3.06 


-2.88 


1.8 


-1.03 


-1.12 


BS 16467-062 


5200 


2.50 


3.07 


-3.77 


-3.58 


1.6 


-1.96 


-1.99 


BS 16477-003 


4900 


1.70 


2.29 


-3.36 


-3.16 


1.8 


0.06 


-0.19 


BS 17569-049 


4700 


1.20 


1.74 


-2.88 


-2.70 


1.9 


0.17 


-0.18 


CD-38:245 


4800 


1.50 


2.33 


-4.19 


-3.91 


2.2 


-0.72 


-0.82 


CS22 169-035 


4700 


1.20 


1.79 


-3.04 


-2.84 


2.2 


-0.33 


-0.30 


CS22 172-002 


4800 


1.30 


2.14 


-3.86 


-3.58 


2.2 


-1.40 


-1.60 


CS22 186-025 


4900 


1.50 


2.10 


-3.00 


-2.80 


2.0 


0.61 


0.41 


CS22 189-009 


4900 


1.70 


2.33 


-3.49 


-3.28 


1.9 


-1.01 


-1.07 


CS22873-055 


4550 


1.00 


1.48 


-2.99 


-2.83 


2.2 


-0.16 


-0.30 


CS22873-166 


4550 


1.00 


1.47 


-2.97 


-2.81 


2.1 


0.07 


-0.07 


CS22878-101 


4800 


1.30 


1.98 


-3.25 


-3.02 


2.0 


-0.41 


-0.36 


CS22885-096 


5050 


2.60 


2.99 


-3.78 


-3.65 


1.8 


-1.44 


-1.47 


CS22891-209 


4700 


1.00 


1.71 


-3.29 


-3.05 


2.1 


0.15 


-0.01 


CS22892-052 


4850 


1.60 


2.14 


-3.03 


-2.85 


1.9 


0.41 


0.20 


CS22896-154 


5250 


2.70 


3.03 


-2.69 


-2.58 


1.2 


0.32 


0.18 


CS22897-008 


4900 


1.70 


2.30 


-3.41 


-3.21 


2.0 


-0.39 


-1.20 


CS22948-066 


5100 


1.80 


2.46 


-3.14 


-2.92 


2.0 


0.46 


0.25 


CS22952-015 


4800 


1.30 


2.02 


-3.43 


-3.19 


2.1 


-0.96 


-1.07 


CS22953-003 


5100 


2.30 


2.68 


-2.84 


-2.71 


1.7 


0.13 


-0.01 


CS22956-050 


4900 


1.70 


2.29 


-3.33 


-3.13 


1.8 


-0.47 


-0.42 


CS22966-057 


5300 


2.20 


2.72 


-2.62 


-2.45 


1.4 


-0.33 


-0.44 


CS22968-014 


4850 


1.70 


2.31 


-3.56 


-3.36 


1.9 


-1.81 


-1.69 


CS2949 1-053 


4700 


1.30 


1.85 


-3.04 


-2.86 


2.0 


-0.15 


-0.36 


CS29495-041 


4800 


1.50 


1.98 


-2.82 


-2.66 


1.8 


-0.23 


-0.36 


CS29502-042 


5100 


2.50 


2.87 


-3.19 


-3.07 


1.5 


-1.99 


-2.05 


CS295 16-024 


4650 


1.20 


1.74 


-3.06 


-2.88 


1.7 


-0.47 


-0.61 


CS295 18-051 


5200 


2.60 


2.93 


-2.69 


-2.58 


1.4 


0.06 


-0.07 


CS30325-094 


4950 


2.00 


2.47 


-3.30 


-3.14 


1.5 


-2.38 


-2.47 


HD2796 


4950 


1.50 


1.95 


-2.47 


-2.32 


2.1 


-0.23 


-0.24 


HD 186478 


4700 


1.30 


1.74 


-2.59 


-2.44 


2.0 


0.05 


-0.12 



Notes. ( ' The synthesised LTE abundances can be found in the First Stars papers; Dwarfs from Bonifacio et al. (2009, B09), and giants from 
Francois et al. (2007, F07) 



10. As discussed in Bergemann et al. (2012a), the model atom of NLTE abundance corrections. For the resonance Srn line at 
by Andrievsky et al. (201 1) is less complete than ours, and we 4077 A, our NLTE corrections are mildly negative for any log g 
incorporated new atomic data, which influences the magnitude and T e /f at [Fe/H] = -3, whereas Andrievsky et al. (2011) ob- 
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Fig. 9: Upper figure: LTE Sr synthesised (Syn) abundances for 
F07, B09 (open circles), and our sample (filled circles - giants 
big, blue; dwarfs small, red). The figure in the middle shows 
the EW-based abundances for the same stars. Lower figure: The 
same samples but with NLTE corrected stellar parameters, and 
NLTE corrected EW-based Sr abundances. 



tain large positive corrections for dwarfs and negative ANLTE 
for giants. 

5. Discussion — Chemical evolution of Sr 

Since we wish to assess the impact of the LTE assumption vs 
NLTE on the chemical evolution of Sr, we selected a handful 
of Sr yields, covering both s- and r-process contributions. We 
probed how LTE vs NLTE abundances behave in a Galactic 
chemical evolution scheme. The yields will briefly be outlined 
below. 

5. 1 . Theoretical predictions of stellar Sr yields 

Here we consider the weak r-process yields from Arcones & 
Montes (201 1) and Wanajo et al. (201 1), and the s-process yields 
from Bisterzo et al. (2010), Frischknecht et al. (2012), and our 
own AGB yields based on the calculations presented in Karakas 



«L -1 




-3 -2 
[Fe/H] 



Fig. 10: Difference in [Sr/Fe] ratios between NLTE (blue dia- 
monds) and LTE (open circles). 



et al. (2012) and Lugaro et al. (2012) that extend down to low 
metallicities of [Fe/H] = -2.3. 

The yields from Wanajo et al. (201 1) describe the ejecta from 
low-mass (~ 9M Q ) faint core-collapse electron-capture super- 
novae (ECSN). These may occur frequently even at low metal- 
licity (e.g. Langer 2012), and their yields can therefore not be 
neglected when considering the evolution of Sr. The amount of 
Sr injected into the interstellar medium (ISM) is, based on the 
self-consistent two-dimensional ECSN model of Wanajo et al. 
(201 1), 1.79 -10~ 4 Mo, which is obtained from the ejecta with an 
entropy ranging from 10 to 25kfi/baryon and an electron frac- 
tion (Y e ) between 0.4 and 0.56. This corresponds to the standard 
yield (W sla „d)- Low Y e values are expected in low-mass progeni- 
tors due to their fast explosions. Convective bubbles expand fast 
enough to inhibit the neutrino reaction that increases the electron 
fraction. However, owing to the self-consistency of the explosion 
model there are no free parameters in the simulations once the 
progenitor model and physics input are chosen. The uncertain- 
ties therefore stem from the systematic uncertainties associated 
with the input to the SN model and possible resolution limita- 
tions, but mainly from the progenitor mass range we adopt in 
the GCE models, which is what gives rise to Wlow and Whigh. 
Here we have selected two mass intervals for Wlow and Whigh 
that are representative of the SN mass used for these SN model 
calculations (see also Langer (2012); Nomoto (1987)). 

Neutrino-driven winds following immediately after super- 
nova explosions, also from more massive progenitors, will like- 
wise contribute to the amount of Sr in the ISM. The wind pre- 
dictions we have incorporated here are based on the compu- 
tations presented in Arcones & Montes (2011). Here we have 
tested the impact that wind parameters, such as entropy, elec- 
tron fraction, and expansion time scale, have on the neutron- 
capture nucleosynthesis in the wind. The effect of progenitor 
mass and progenitor metallicity remains an open question, since 
the supernova models still have uncertainties that are too large to 
constrain these quantities. Nonetheless, we note that since this 
process is a primary process, the impact of metallicity is not 
the most important factor when trying to constrain the yields. 
Assuming a neutron-rich wind with the electron fraction con- 
stricted to 0.4 < Y e < 0.49 (as currently suggested by Martinez - 
Pinedo et al. 2012; Roberts & Reddy 2012; Roberts 2012), we 
can try to loosely confine some parameters to realistic ranges by 
setting an entropy interval of 50 to 150 £#/baryon and a wind ex- 
pansion time scale limited to a few milliseconds. These are typ- 
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ical values found in hydrodynamical wind simulations (Arcones 
et al. 2007; Fischer et al. 2010) and lead to Sr yields spanning 
10~ 4 - 1O~ 7 M , where the largest contribution comes from 12 
- 25 M supernovae, and the smallest yield could be assigned 
to 8 - 12 M SN explosions. Generally speaking, the smaller 
Sr yields can be assigned to low-mass progenitors, while more 
massive ones produce larger Sr yields. 

The first s-process yield we consider here are those of 
Bisterzo et al. (2010), who provide the yields from AGB stars 
in the mass range 1.3 -2M at different metallicities ([Fe/H]): 0, 
-0.8, -1.6, and -2.6 (S. Bisterzo, priv. comm.). These yields have 
been calculated with the FRANEC (Frascati Phapson-Newton 
Evolutionary Code) that uses reaction rates from the KADoNiS 
and NACRE databases. The neutron-capture elements are cre- 
ated in 13 C pockets and brought to the stellar surface during 
thermal pulses. The AGB star experiences a mass loss in the 
range 10~ 4 - 10~ 7 M G /yr where the interval between successive 
thermal pulses lasts for about 10 4 — 10 5 years depending on the 
AGB core mass (e.g., see model data published in Karakas et al. 
2012 or Cristallo et al. 2011). The impact at low metallicity 
([Fe/H]< -2.6), which we are interested in, is minor and would 
need to be extrapolated. However, yields from low-mass AGB 
stars are important at higher metallicities. 

We also use AGB yields based on the nucleosynthesis cal- 
culations presented in Karakas et al. (2012) and Lugaro et al. 
(2012). These calculations use reaction rates from the NACRE 
and JINA databases, which includes the KADoNIS neutron- 
capture cross sections. In comparison to the yields from Bisterzo 
et al. 2010, these span a range from 0.9 M to 6M at [Fe/H] = 
-2.3; at higher metallicities we cover a mass range from 1.25M 
to 6M at [Fe/H] = -0.15 and 5M to 8M at [Fe/H] = +0.14 
(slightly super solar). For the heaviest of these stars, the main 
neutron source is the 22 Ne(a,n) 25 Mg reaction, and a 13 C pocket 
has been introduced. However, for the lowest mass stars, we in- 
troduce some partial mixing of protons at the deepest extent of 
each third dredge-up episode in order to form a 13 C pocket (see 
detailed discussion in Lugaro et al. 2012). For most models the 
mass-loss rate used on the AGB is the semi-empirically derived 
Vassiliadis & Wood (1993) mass loss formulae. One of the ma- 
jor uncertainties in AGB modelling is the mass-loss rate, espe- 
cially at the lowest metallicities, and for this reason we exper- 
iment with variations to the mass loss in the intermediate-mass 
AGB models. We refer to the discussions in Karakas (2010) and 
Karakas et al. (2012) for more details. 

The s-process yields of massive rapid rotating stars are from 
Frischknecht et al. (2012); these results are more recent than 
Pignatari et al. (2008), and Frischknecht et al. (2012) based their 
yield predictions on a stellar evolution code. The code used is 
the Geneva stellar evolution code (GENEC), with reaction li- 
braries from REACLIB, and neutron captures from KADoNiS. 
Moreover, the results by Frischknecht et al. (2012) are based 
on the still favoured reaction rate for ll O(a, y) by Caughlan & 
Fowler (1988). We explore the uncertainty linked to this reac- 
tion rate decreasing it by a factor of 10. In this way we can 
generate upper and lower limits for what we might expect as 
Sr yields from these massive rotating stars. In general the above- 
mentioned s-process yields are in the range 10~ 6 - 1O~ 9 M . (For 
comparison the non-rotating stars barely produce any Sr.) 

5.2. The impact of different stellar yields on Galactic 
chemical evolution (GCE) models 

We now show a chemical evolution model computed with the 
different assumptions for the stellar yields discussed above. 



Because the goal of this section is to check for the broad im- 
pact of the different stellar yields on the chemical enrichment, 
we adopt a homogeneous chemical evolution model to guide our 
discussion. For that we use the halo chemical evolution model of 
Chiappini et al. (2008), which we briefly summarise here: 

1) the infall law of the primordial gas follows a Gaussian func- 
tion, 

2) the occurring outflow from the system is proportional to the 
star formation rate (SFR). 

The timescale for the formation of the halo is fast (less than 
0.5 Gyr). The other crucial element for our study is iron, for 
which we adopt the predictions from Woosley & Weaver 1995, 
(WW95). As shown in Cescutti & Chiappini (2010), this com- 
bination of parameters is able to produce a synthetic metallicity 
distribution function, which is in good agreement with the one 
observed in the Galactic halo (Lai et al. 2008; Schorck et al. 
2009). As a result we assume that this model follows the correct 
timescale for the chemical enrichment. 

We have not run individual models for each set of AGB nu- 
cleosynthesis prediction. The reason is that, even though our 
chemical evolution model calculates the enrichment from low- 
mass AGB stars and SNe la, only the more massive stars (SN II) 
can be considered drivers of the early chemical enrichment in the 
halo. This is due to the considerable longer timescales needed to 
evolve AGBs and SNe la. We have only tested the metal-poor 
yields for Sr (i.e. from the metal-poor extension to the yields 
from Lugaro et al. 2012), for which we have predictions up to 
6M stars (compared to the 2 - 3M from Bisterzo et al. 2010; 
Cristallo et al. 201 1). The massive AGB stars are more likely to 
contribute to the chemical enrichment during the formation of 
the halo, i.e. at lower metallicities. Nevertheless, the enrichment 
produced by the 6M AGB stars do not influence the overall re- 
sults. Only when extremely low yields from even more massive 
stars (SN II) are adopted (lower than Wlow), is it possible to see 
the influence of the heavy AGB yields. 

We have computed different chemical evolution models with 
different nucleosynthesis assumptions, namely: 

- Case 1) Wanajo low (Mow): 1.79-1O~ 4 M in the range 
9.5-lOM ; 

- Case 2) Wanajo high (Whigh): 1.79-1O~ 4 M in the range 
8-10M o ; 

- Case 3) Arcones & Montes low (AMlow): 1O~ 7 M mass 
range 8-25M ; 

- Case 4) Arcones & Montes high (AMhigh): 1O~ 4 M in mass 
range 8-25M 3 ; 

- Case 5) Frischknecht low (Flow) with rotation and their 
'standard' value for the reaction rate, generalising the pro- 
duction of the 25M in their paper to a range of 15-4OM ; 

- Case 6) Frischknecht high (Fhigli) calculated as Flow but 
with a decreased value for the reaction 17 rate (producing 
higher Sr yields) 

In Fig. 11a) (top) we show the models Wlow and Whigh, 
which differ in the level of Sr enrichment by the ECSN due to 
the supernova mass used in the GCE model. The high level of 
Sr production in a narrow range of masses produces in the case 

3 Arcones & Montes standard (Astand): 1O~ 4 M in mass range 12- 
25M e and 1CT 7 M G in mass range 8-12M G and M> 25M e . 
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Fig. 11: GCE model predictions compared to Case 1+2), which 
is based on ECSN yields from S. Wanajo in the top panel 
a). The predictions are in the middle panel b) compared to 
Case 3+4) which is neutrino-driven winds from Arcones & 
Montes, and in the bottom panel c) the comparison is made 
to Case 5+6), namely fast-rotating stars from U. Frischknecht. 
The dashed/solid lines represent the upper/lower limits to the Sr 



of Whigh an important bump in the [Sr/Fe] predictions. This 
roughly appears as a knee in the [Sr/Fe] trend at [Fe/H]~ -3. 
However, the Whigh model produces too much Sr compared to 
the observed values. Including CEMP-s stars, we would have 
found some stars closer to the Whigh model only at the lowest 
metallicity. The Wlow model predicts [Sr/Fe] ratios very close to 
the observations at all [Fe/H] values. 

Model AMIow and AMhigh are presented in Fig. 1 lb) (mid- 
dle). In these models the Sr is produced in a wide range of mas- 
sive stars, which leads to a rather flat [Sr/Fe] trend. The nearly 3 
dex difference between the upper and lower curves {AMIow and 
AMhigh) reflects the variations and uncertainties in the Sr yields. 

From the results obtained from the two pairs of models we 
cannot draw strong conclusions as to the nucleosynthetic origin 
of Sr in our metal-poor star sample. The only message we can 
read from them is that some combination of two r-process sites, 
can roughly explain the observationally derived abundances. At 
this level, the observations/data can still be used to constrain the 
theoretical yields, rather than the other way around. 

Finally, from the two models generated using the 
Frischknecht et al. (2012) yields (Fig. 11c) bottom panel), we 
can conclude that the s-process from fast-rotating massive stars 
may be an important source of Sr during the Galactic halo forma- 
tion (see Cescutti et al. 2013). The early s-process taking place in 
metal-poor, fast rotating, massive stars should be coupled/added 
to at least one r-process. 

Since the difference between the LTE and NLTE Sr abun- 
dances is small, despite the slight horizontal shift caused by 
the positive NLTE effect on metallicity, the methodology does 
not affect the interpretations in a GCE context (see Fig. 12). 
Therefore only the same GCE conclusions can be drawn under 
LTE and NLTE. We stress that this is only true for an element 
like Sr. Owing to the combined NLTE effects in Sr and Fe, there 
are no stars below [Fe/H] =-3.5 with high strontium abundances 
([Sr/Fe]> 0). Our NLTE corrected abundances still show a large 
star-to-star scatter as is also found in Andrievsky et al. (201 1). 

Without precise theoretical yields it is impossible to draw 
firm conclusions based on the abundance measurements. On the 
other hand, the presence of the spread can be explained in a 
stochastic chemical evolution scenario. 

Many advances have been made in modelling AGB stars, 
SN, and the r-process, as well as significant improvements in 
the determination of atomic data for heavy elements. These have 
in turn improved the GCE models. However, the r-process yields 
still need to be better constrained not to span 3 dex. Until these 
improvements are made, yield predictions should only be used as 
guiding upper and lower limits. The yields from the fast-rotating 
stars and the AGB stars also face challenges when trying to con- 
strict the 13 C pockets, pulse duration, mass loss, and poisons in 
order to improve the networks and, in turn, yields. Only within 
the last decade have SN models managed to explode, and the 
treatment of 13 C pockets and mass loss from AGB stars is im- 
proving (e.g. Matsuura et al. 2007; Mattsson et al. 2008; Sloan 
et al. 2012). For example, by using AGB stars in star clusters 
it is possible to constrain uncertain parameters related to con- 
vection and mass loss (e.g. Lebzelter et al. 2008; Kamath et al. 
2012), while 13 C pocket sizes can be constrained using data from 
a variety of observational data (e.g., Lugaro et al. 2003; Bonacic 
Marinovic et al. 2007; Bisterzo et al. 2012). 

There are still many challenges, including a detailed under- 
standing of the formation mechanism of 13 C pockets, which is 
currently unknown, and a significant uncertainty in s-process 
models of AGB stars (e.g., see discussion in Herwig 2005). 
The yields from fast-rotating massive stars also face some of 
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Fig. 12: Comparison of the best homogeneous GCE model (solid 
line) to Sr LTE abundances (open circles) in the upper panel a), 
and to NLTE Sr abundances (filled blue diamonds) in the lower 
panel b). 



the same challenges including an accurate description of mass 
loss and convection in stellar models (e.g., see review by Langer 
2012). 

However, on the basis of Figs. 1 1 and 12 we believe that at 
least two different sources (e.g. neutrino-driven winds from mas- 
sive stars or low-mass ECSN) are needed to explain the [Sr/Fe] 
abundances we derived under both LTE and NLTE. An early 
s-process could take place in fast-rotating stars and provide a 
possible explanation for the low [Sr/Fe] we obtain at very low 
metallicity. 



Cross-disciplinary work is crucial for computing better yields 
and, in turn, better chemical evolution models. Only then will 
such calculations allow us to accurately predict the behaviour 
and evolution of heavy elements such as Sr. This information is 
needed in the era of high-resolution surveys, where we have a 
large flow of data coming in. We are now capable of analysing 
the stars and extracting very accurate abundances with uncer- 
tainties < 0.25 dex, compared to what can be obtained from nu- 
cleosynthetic yields and evolved stellar models. This indicates 
that observationally derived abundances are more likely to con- 
strain the parameter space of the yield predictions and the GCE 
models than vice versa. 

On the basis of this study we can conclude the following 
about Sr in metal-poor stars. 

- When deriving Sr I abundances, NLTE corrections should al- 
ways be applied in order to obtain ionisation equilibrium. 

- In the interval -3 .0 < [Fe/H] < 1.0 the chemical evolution 

of Sr, as derived from the SrII resonance lines, is similar 
under LTE and NLTE. However, below [Fe/H] = -3.0, the 
NLTE corrections to Srn lines are important to obtain the 
correct Sr abundances. 

- The abundances obtained from the Sr n lines are sensitive to 
surface gravity. If the latter parameter is obtained from the 
ionisation equilibrium of Fe, LTE approximation should not 
be used, because it leads to large systematic errors in log g 
of up to +0.8 dex. 

- In their current state the Sr yields are too uncertain to clearly 
disentangle contributions from different processes or sites. 
We may instead use the observationally derived abundances 
to constrain the parameter space of the model predictions. 

In summary, it is not sufficient to account for NLTE ef- 
fects in the line formation of the Sr lines. NLTE effects must be 
taken into account in determining stellar parameters, i.e., surface 
gravities and metallicities. Alternatively, abundances determined 
from the lines of the majority species, Sr n and Fe n, can be used. 
The LTE assumption is a trustworthy chemical evolution tracer, 
in the interval -3 < [Fe/H] < 0, for the Sr II abundances from 
dwarfs calculated with gravities stemming from parallaxes and 
temperatures based on accurate photometry. The metal-poor gi- 
ants, which are the best targets at low metallicity, are biased 
by the LTE assumption. The parameters of the giants and their 
abundances must be computed under NLTE. 

With the current uncertainties on the stellar yields, which 
span two to three orders of magnitude, we cannot draw strong 
conclusions on the chemical evolution of Sr in the early Galaxy. 
Neither can we precisely extract the various sites that contribute 
to the creation of the large star-to-star scatter, although at least 
two sites seem necessary. The yield predictions and the GCE 
model used in this work provide both upper and lower limits to 
the highly scattered stellar strontium abundances. 
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-0.02 


0.12 


< -0.35 


0.28 


< -0.07 


HD 19445 


5982. 


4.38 


-2.13 


-0.05 


0.13 


< 0.14 


4.38 


-2.10 


-0.12 


0.11 


-0.05 


0.06 


< 0.05* 


0.34 


< 0.39* 


HD 10603 8 


5950. 


4.33 


-1.48 


0.65 


0.5 


0.33 


4.33 


-1.45 


0.6 


0.5 


-0.02 


0.48 


0.3 


0.3 


0.6 


HD121004 


5711. 


4.46 


-0.73 


0.21 


0.15 


0.2 


4.46 


-0.71 


0.2 


0.15 


-0.01 


0.14 


0.14 


0.24 


0.38 


HD122196 


6048. 


3.89 


-1.81 


0.05 


0.24 




3.89 


-1.75 


-0.01 


0.22 


-0.03 


0.19 








HD 122563 


4665 


1.65 


-2.60 


-0.09 


-0.05 


< -0.6 


1.65 


-2.50 


-0.09 


-0.05 


0.0 


-0.05 


< -0.7 


0.45 


< -0.25 


HD 140283 


5777 


3.70 


-2.58 


-0.33 


-0.15 




3.70 


-2.38 


-0.54 


-0.36 


-0.01 


-0.37 








HD175305 


5100. 


2.70 


-1.38 


-0.02 


0.1 


-0.35 


2.70 


-1.34 


-0.06 


0.11 


0.0 


0.11 


-0.4 


0.37 


-0.03 


G6412 


6464 


4.30 


-3.24 


-0.05 


0.0 




4.30 


-3.12 


0.04 


0.05 


0.12 


0.17 








G6437 


6494. 


3.82 


-3.17 


-0.09 


0.08 




4.23 


-3.00 


-0.05 


0.08 


0.09 


0.17 








BD133442 


6450. 


4.20 


-2.56 


0.26 


0.3 




4.42 


-2.47 


0.14 


0.27 


-0.06 


0.21 








CS303 12-059 


5021. 


1.90 


-3.06 


0.46 


0.5 




2.41 


-2.89 


0.28 


0.35 


-0.04 


0.31 








CS31082-001 


4925. 


1.51 


-2.81 


0.59 


0.7 




2.05 


-2.63 


0.65 


0.62 


-0.02 


0.6 








HD74462 


4590. 


1.84 


-1.48 


0.1 


0.0 


-0.5 


1.98 


-1.43 


0.04 


-0.12 


0.0 


-0.12 


-0.5 


0.34 


-0.16 


HD 12623 8 


4900 


1.80 


-1.92 


0.02 


0.0 


-0.34 


2.02 


-1.85 


0.03 


0.05 


0.0 


0.05 


-0.46 


0.43 


-0.03 


HD126587 


4950. 


1.90 


-3.01 


0.39 


0.37 


< -0.05 


2.36 


-2.86 


0.34 


0.25 


-0.04 


0.21 


< -0.15 


0.4 


< 0.25 


HE03 15+0000 


5050. 


2.05 


-2.81 


0.17 


0.39 




2.47 


-2.67 


0.03 


0.26 


-0.03 


0.23 








HE1219-0312 


5100. 


2.05 


-2.99 


0.22 


0.29 




2.58 


-2.81 


0.08 


0.17 


-0.05 


0.12 
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Notes. w Larger uncertainty ±0.1 dex in measurement 
i<} upper limit on the abundance 



